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It is commonly believed that our Universe has experienced at least two different stages of acceler- 
ated expansion. The early stage is known as inflation and the current acceleration is driven by dark 
energy. Observing inflation and dark energy dynamics are among the most important aspects of the 
current cosmological research. In this paper we make a first detailed probe of the possible degenera- 
cies between dynamical inflation and dark energy in light of the current cosmological observations. 
We have combined Type-la supernova (Riess "gold" samples), galaxy clustering (SDSS 3D power 
1 spectra) and the Wilkinson Microwave Anisotropy Probe (WMAP) observations and performed a 

q ' global analysis using the Markov chain Monte Carlo (MCMC) method. We find the inclusion of 

, inflation and dark energy parameters together make the parameter spaces broader and the degen- 

eracy between the inflation dark energy parameters are nontrivial: the allowed/preferred behavior 
of dynamics in one sector is dependent on the prior of the other sector. Interestingly a devia- 
tion from scale-invariant primordial spectrum is sightly more preferred by the current cosmological 
observations when one marginalizes over dynamical dark energy models. 
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I. INTRODUCTION 
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ff) i Recent advances in both theoretical and observational cosmology have revealed that our Universe has experienced 
two different stages of accelerated expansion. One is the inflation in the very early universe when its tiny patch 
was superluminally stretched to become our observable Universe today P, Q|. It can naturally explain why the 
universe is flat, homogeneous and isotropic. Inflation is driven by a potential energy of a scalar field called inflaton 
and its quantum fluctuations turn out to be the primordial density fluctuations which seed the observed large-scale 
structures (LSS) and anisotropy of cosmic microwave background radiation (CMB). In the past decade, inflation 
theory has successfully passed several nontrivial tests. In particular, the released first year Wilkinson Microwave 



Anisotropy Probe (WMAP) data |3j have detected a large-angle anti-correlation in the temperature-polarization 
i cross-power spectrum, which is the signature of adiabatic superhorizon fluctuations at the time of decoupling [ij. The 
cosmological observations such as the first year WMAP(WMAPl) and the Sloan Digital Sky Survey (SDSS) are in 
good agreement with adiabatic and scale invariant primordial spectrum(see e.g. which is consistent with single 

field slow rolling inflation predictions. On the other hand, however, inflation generically predicts primordial spectrum 
with some deviations from scale invariance and it is crucial to probe the scale dependence of primordial spectrum in 
light of the cosmological observations. The tentatively interesting features of the first year WMAP observations have 
aroused a lot of interests to explain such effects with the dynamics of inflation 0. 

In 1998 the analysis of the redshift-distance relation of Type la supernova (SNIa) has revealed the existence of 
the second stage of accelerated expansion that has started rather recently when a mysterious new energy component 
dubbed dark energy (DE) has dominated the energy density of the Universe. 0,0. Recent observations of SNIa have 
confirmed the accelerated expansion at high confidence level |ToL fill fl2L fl3| . The nature of dark energy is among the 
biggest problems in modern physics and has been studied widely. A cosmological constant, the simplest candidate 
of DE where the equation of state (EOS) w remains -1, suffers from the well-known fine-tuning and coincidence 
problems. Alternatively, dynamical dark energy models with the rolling scalar fields have been proposed, such as 
quintessence [hi Il5j . phantom 0] and k-essence[3>[^. Given that currently we know very little on the theoretical 
aspects of dark energy, the cosmological observations play a crucial role in our understanding of dark energy. If the 
observations strongly favor the dark energy component with a deviation from w = — 1, then the cosmological constant 
which puzzles the theorists would not be the source driving current acceleration of the universe. 

Since there are many common aspects of the background physics, albeit the striking difference of the energy scale 
between the two stages of accelerated expansion, it would be nice if one could unify inflation with dark energy. 
The model of quintessential inflation[19j has been proposed in history trying to unify the two epoches of accelerated 
expansions. The first year WMAP observations show the lack of the temperature-temperature correlation power on 
the largest scales and many studies have attributed it to the suppressed primordial spectrum due to some dynamics 
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of inflaton. Intriguingly, in Ref.|2(| the authors have attributed the lack of power to the isocurvature fluctuations in 
the quintessence which may be generated during inflation, thus the different dynamics of the dark energy sector and 
inflation can lead to similar effects on the observations, where such a degeneracy may reflect the possible connections 
between dark energy and inflation. 

Observing inflation and dark energy dynamics are among the most important aspects of the current cosmological 
research. So far to our knowledge in the previous observational studies of dynamical dark energy or inflation not 
enough attention has been paid to the possible correlation between the two sectors: dark energy and inflation. In this 
paper we make a first detailed probe of the possible degeneracies between dynamical inflation and dark energy in light 
of the current cosmological observations. Our paper is structured as follows: in Section II we describe the method 
and the data; in Section III we present our results by global fittings; finally we give our summary and conclusions in 
Section IV. 



II. METHOD AND DATA 



The role of the primordial scalar spectral index n s in the probe of inflation is somewhat similar to that of dark 
equation of state w in the measurement of dark energy. The case of n s = 1 corresponds to the scale invariant 
Harrison-Zel'dovich spectrum, where the scalar spectrum is constant over all scales and w = — 1 corresponds to the 
cosmological constant where the density of dark energy does not evolve with time. A detection of deviation from 
n s = 1 or w = — 1 would mark a breakthrough in our understanding of inflation and dark energy. To make our study 
on inflation and dark energy symmetric and in some sense model independent, we parametrize n s and w with linear 
expansions and keep them to the first order: 

n s (k) = n s (k s0 ) + a s In > (!) 

where k S Q is a pivot scale which is arbitrary in principle and a s is a constant characterizing the "running" of the scalar 
spectral index. The scalar spectral index n s is related to the primordial scalar power spectrum V x (k) by definition: 

n s (k) = dV x {k)/d\nk + 1. (2) 

Correspondingly V x (k) is now parametrized as[2lj: 

lnV x (k) = lnA s + (n s (k s0 )-l)ln(J^) +^ ( ln (7^)) ■ ( 3 ) 
Accordingly for the dark energy sector the equation of state is parametrized as |22j: 

wde{o) — w + wi(l - a) (4) 

where a = 1/(1 + z) is the scale factor and w\ characterizes the "running" of the equation of state. The equation of 
state is defined as the pressure over the energy density, which leads to the evolution of dark energy density via the 
energy conservation of dark energy: 

PQ (a)/ PQ (a - 1) = a -3(i+-o+- 1 ) exp [3wi{a _ x)] _ (5) 

The current publicly available codes like CMBfast[13, and CAMBjii, have well adapted the constant 
running of the prim ordial spectral index. For the constant running of the dark energy equation of state it is not so 
straightforward 27, 28] . For the common dynamical dark energy models like quintessence or phantom, the equation 
of state does not get across -1, and for the global fittings through parametrizations of w the dynamics of quintessence 
or phantom has been fully enclosed in w in which sense one can even reconstruct the potentials of quintessence 
and phantom with w. However in the cases where w gets across -1 during evolution, single field scalar dark energy 
models like quintessence, phantom or k-essence cannot realize such a transition and on the other hand, for the 
parametrizatized dark energy one is encountered with the problem of divergence for the perturbation equations |28|. 
These properties have been investigated in details in our companion paper in Ref.|28|. Dark energy with equation of 
state getting across -1 has been dubbed quintom in the sense that it resembles the combined behavior of quintessence 
and phantom HI EH We have shown that in general to realize the model of quintom one needs to add extra 
degrees of freedom beyond the standard single field case[2£|. In the extant two-field-quintom model and the single 



3 



field model with a high derivative term|32(, the perturbation of dark energy is shown to be continuous when w gets 
across -1. 

For the parametrization of the equation of state which gets across -1, we introduce a small positive constant e to 
divide the full range of the allowed value of w into three parts: 1) w > — 1 + e; 2) — 1 + e > w > — 1 — e; and 3) 
w < — 1 — e. Working in the conformal Newtonian gauge, one can easily describe the perturbations of dark energy as 
follows |33] ]: 

5 = -(l + w)(6-3<f>)-m(c 2 s -w)5 , (6) 

9 = -H(l-3w)0 —8 + k 2 (^- + V) . (7) 

1 + w 1 + w 

Neglecting the entropy perturbation contributions, for the regions 1) and 3) the equation of state does not get 
across —1 and perturbations are well defined by solving Eas. i|6l7[l . For the case 2), the perturbation of energy density 
5 and divergence of velocity, 9, and the derivatives of 5 and 9 are finite and continuous for the realistic quintom dark 
energy models. However for the perturbations of the parametrized quintom there is clearly a divergence. In our study 
for such a regime, we match the perturbation in region 2) to the regions 1) and 3) at the boundary and set|28| 

6 = , = 0. (8) 

In our numerical calculations we've limited the range to be \Aw = e\ < 10~ 5 and we find our method is a very 
good approximation to the multi-field quintom. For more details of this method we refer the readers to our previous 
companion paper |2j| ■ 

In this study we have implemented the publicly available Markov Chain Monte Carlo package CosmoMCjlil], which 
has been modified to allow for the inclusion of Dark Energy perturbation with the equation of state getting across 
— 1. We assume purely adiabatic initial conditions and a flat Universe. Our most general parameter space is: 

P = (u b ,uj Cl <d s ,T, w ,Wi,n a ,a s ,ln(10 w A s ),r) (9) 

where uif, = fli,h 2 and ui c = Vl c h 2 are the physical baryon and Cold Dark Matter densities relative to the critical 
density, S is the ratio (multiplied by 100) of the sound horizon to the angular diameter distance at decoupling, r is 
the optical depth to re- ionization, A s , n s and a s characterize the primordial scalar power spectrum, r is the tensor to 
scalar ratio of the primordial power spectrum and for the tensor slope tit we have fixed r = Snx- For the pivot of 
the primordial spectrum we set k s0 — 0.05Mpc _1 . Furthermore, we make use of the Hubble Space Telescope (HST) 
measurement of the Hubble parameter Hq = lOOh km s _1 Mpc^ SS] by multiplying the likelihood by a Gaussian 
likelihood function centered around h = 0.72 and with a standard deviation a = 0.08. We also impose a weak Gaussian 
prior on the baryon and density £!/,/i 2 = 0.022 ± 0.002 (1 a) from Big Bang Nucleosynthesis [3fij . 

In our calculations we have taken the total likelihood to be the products of the separate likelihoods (£,) of CMB, 
LSS and SNIa. In other words defining \f = — 2 log A, we get 

Xtotal = XCMB + XLSS + XsNIa ■ (10) 

We have included the first-year temperature and polarization data with the routine for computing the likelihood 
supplied by the WMAP teamj3,E3. We use the code developed in Ref . [38j| to fit the 3D power spectrum of galaxies 
from the SDSS and the bias factor of SDSS has been used as a continuous parameter to give the minimum \ 2 value. 
For the SNIa data we use the 157 gold-sample data published by Riess et a/|ll| and in the calculation of the likelihood 
we have marginalized over the nuisance parameter 39]. 

Compared with Riess "gold" sample of SNIa, the first year Supernova Legacy Survey (SNLS) sample^] has 
been measured with one telescope and the systematics are better constrained. However for the time being their 
determination on cosmological parameters are marginal comparable |40j . Due to the fact that the two samples of SNIa 
have been from different methods of light curve fitting a combination of them is not so straightforward. For LSS for 
the time being the Two degree Field (2dF) 01 an d SDSS [3 are also comparable and for the current study we have 
not combined the two data sets simultaneously. 

For each regular calculation, we run 8 independent chains comprising of 150, 000 — 300, 000 chain elements and 
spend thousands of CPU hours to calculate on a supercomputer. The average acceptance rate is about 40%. We 
test the convergence of the chains by Gelman and Rubin criteria an d find R — 1 is of order 0.01 which is more 
conservative than the recommended value R—l < 0.1. 
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TABLE 1. Median values and la constrains on n s , r, a s , wo and wi for the three models as explained in the text, shown 
together with the minimum \ 2 values . For r the 2a upper bounds have been shown instead. 
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III. RESULTS 

To get the clear effects on the possible correlations between dynamical inflation and dark energy we mainly focus 
on the three kinds of models: one is the ACDM cosmology where n s , r and a s are enclosed in the inflationary sectors, 
one is with scale-invariant primordial spectrum (r not included) where wq and w\ are included in the dark energy 
sector, and the third one includes all of the free parameters like n s , r, a s , wq and wi . 

Table 1 lists all of the relevant 1-dimensional median values and la constrains for the three models specified above. 
Shown together are the minimum \ 2 values for each case. For the constraints on r only 2a upper bounds have been 
shown. The effects of different priors on the inflationary and dark energy sectors will be explained in details with the 
2-dimensional graphics of the sub parameter space below. 




FIG. 1: Constraints on the primordial spectrum parameters by marginalizing over ACDM cosmology (Blue Dash line) and 
over dynamical dark energy with equation of state w = too + Wi(l — a) (Red Solid line). The two dimensional contours stand 
for 68% and 95% C.L. 

In Fig. ^we delineate the two dimensional contours of the primordial spectrum parameters by marginalizing over 
ACDM cosmology (Blue Dash lines) and over dynamical dark energy with equation of state in Eq. Q (Red Solid 
lines). The contours stand for 68% and 95% C.L. Although in general the case with dynamical dark energy gives 
less stringent constraints on the primordial spectrum, the ACDM case does not overlap fully with the dynamical 
dark energy case. Intriguingly for the combination with the inclusion of dark energy dynamics, the constraints on 
the tensor contributions exhibit different behaviors at 68% and 95% C.L. And at 68% C.L. the constraint on r is 
more stringent with the inclusion of dark energy dynamics, which is different from naively expected. This can be 
explained by the mild preference on the dark energy dynamics in the combination of WMAPl±SDSS±SNIa(Riess) 



1 For the simplicity of discussions we do not discuss the effects on A s and Qde> although we did vary them in our fittings. 
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and the degeneracy between dark energy dynamics and r on the largest scales of CMB. From Table 1 we can find 
the preference of dynamical DE is slightly larger than la from the different minimum \ 2 values as well as from the 
1-dimensional constraint on the parameter of Wq. Also we can see some difference in the constraints on a s and the 
value of n s at 0.05 Mpc -1 . Note the best fit values of a s in some sense deviate a lot from zero, the difference on the 
probe of scale variance may be non-negligible for the two cases. We get a s — —0.032^® ^® for the ACDM cosmology 
and —0.040^0 036 instead for the case with dynamical dark energy. Roughly speaking a scale invariant primordial 
spectrum is not yet ruled out in the ACDM cosmology while disfavored more than la when we allow the dynamics of 
dark energy! We will turn to this in more details in the later part of this paper. 
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FIG. 2: 68% and 95% constraints on the dynamical dark energy model w = wo + wi(l — a) marginalizing over scale-invariant 
primordial spectrum (Blue Dash line) and over primordial spectrum with a running spectral index (Red Solid line). In the 
latter case tensor contributions are included. The Dotted lines stand for wo = —1 and wo + w i = —1. 



In Fig. [21 we plot the 68% and 95% constraints on the (wq, w\) contour marginalizing over scale- invariant primordial 
spectrum (Blue Dash line) and over primordial spectrum with a running spectral index (Red Solid line). We find 
interestingly the 1 a regions do not overlap fully and the cosmological constant case lies at the edge where the 1 
a regions over lap. Numerically we get wq = —1.28 ± 0.23, w\ = 0.569to'g27 f° r the scale-invariant case and wq = 
— 1.28 ± 0.25, w\ = 0.59llo 881 with a running spectral index. We cut the parameter space of wq — w\ into four areas 
by the line of wq = —1 and wo + w\ = —1, Quintessence, Phantom, Quintom A (wde is Phantom-like today but 
Quintessence-like in the past) and Quintom B (Dark Energy has wde > — 1 today but wde < — 1 at higher rcdshifts). 
Quintom A almost occupies the 1 a region. In a previous companion paper in the probe of dynamical dark energy, 
we got similar results where for the primordial spectrum part we included only n s and j4 s |40|. However inflation 
predicts nonzero tensor contribution and it is a crucial parameter towards verifying the theory of inflation, in the 
current paper we consider r contributions in cases we do not simply assume a Harrison-Zel'dovich spectrum. 

In the probe of deviations from scale invariance and dynamical dark energy it is sometimes more intuitive to get 
constraints on the whole slopes of n s (k) and w(z) on the relevant scales. We have got the constraints on n s (k) and 
w{z) by using the covariance matrix of n s {k s o), a s and wq, w\ from our MCMC results. In Fig. |3|and Fig. 0|we plot 
the constraints on n s (k) and w{z) for the two different cases. In Fig. |3|it is obvious that scale-invariant primordial 
spectrum is more disfavored in the case with dynamical dark energy, although it does allow more parameter space 
compared with the ACDM cosmology. It is intriguing that if we start always from ACDM cosmology in the probe 
of the primordial spectrum information, we might fail to find the true physics which is crucial for understanding 
inflation. Neglecting dark energy dynamics might lead to the bias on the probe of inflation, although the statistical 
significance is rather low. 

In Fig. 0]we find that in both cases dynamical dark energy is favored more than lc compared with the cosmological 
constant. This happens that in both cases wq < —1 is favored more than la, as shown previously. The la regions 
in Fig. ^differ little when marginalizing over scale-invariant primordial spectrum and over primordial spectrum with 
a running spectral index, which can also be easily understood from Fig. [2F and the minimum \ 2 values in Table 1. 



2 Ref.|3 has done some similar analysis on this point, where different observational datasets have been used with somewhat different 
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FIG. 3: 68% and 95% constraints on the slope of n a (k) by marginalizing over ACDM cosmology (Left) and over dynamical 
dark energy with equation of state w = wo + wi(l — a) (Right). The blue thick line stands for the scale-invariant primordial 
spectrum. 
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FIG. 4: 68% and 95% constraints on the evolution of dynamical dark energy model w — wo + wi(l — a) marginalizing over 
scale-invariant primordial spectrum (Left) and over primordial spectrum with a running spectral index (Right). In the latter 
case tensor contributions are included. 



results. They did not include the effects of dark energy perturbations and study the constraints on dynamical dark energy with also a 
and r. In this paper we have tried to probe the deviation from scale variance and dynamical dark energy in equal weights. 
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Our results here are somewhat similar to a previous analysis in Ref. |43j. The mild preference of dynamical dark 
energy over a cosmological constant lies mainly on the Riess sample of SNIa[H. On the other hand such a preference 
also shows imprints on large scales of CMB multipoles via the Integrated Sachs- Wolfe (ISW) effect and dark energy 
perturbations, which are correlated with r and a s . Future improvements on the SNIa observations will also inversely 
help to probe the information of the primordial spectrum, as implied in our Fig. n an d Fig. [21 

IV. SUMMARY AND DISCUSSIONS 

In this paper we have preformed an analysis of global fitting allowing simultaneously the dynamics in both inflation 
and dark energy sector. Our result shows that there is generically some correlation between the inflation and dark 
energy parameters and the parameter space is typically enlarged for the probe of dynamics in either sector compared 
to cases where one assume there can be some dynamics in only one of these sectors. When we allow some dynamics 
in the dark energy sector we find interestingly a deviation from scale-invariant primordial spectrum is sightly more 
favored by current cosmological observations than assuming a ACDM cosmology. 

We should stress that starting from the theoretical aspects it is not so natural to predict exactly a constant running 
on a s or w\, on the other hand given the current precision of the observations and also for a generic study of the 
dynamics we are simply using some of the parametrizations on each sector. In general it is not so easy to achieve 
a precision cosmology |44|, given the fact that for CMB (also LSS and lensing) on large scales we have the cosmic 
variance, on small scales we have the secondary effects which are non-negligible compared with the damped power 
itself and for LSS on small scales we have the not yet well-understood factor of bias and the complicated nonlinear 
evolutions. And for SNIa the calibration is in some sense an open issue. Actually with our better understandings of the 
Universe and the developments of technology we are able to probe the Universe with more versatile observations with 
time going on. Nowadays we are also hoping to extract the information of the Universe through the observations from 
weak and strong lensing, radio galaxies, Gamma-ray bursts(GRB), Lyman-a forest(Lya) and the bayonic oscillations 
inherited in the LSS observations. However we also need to understand better their inherited systematics before 
we can extract more robust and precise information of our Universe. In this paper we have not considered small 
scale CMB observations or those like Lya or GRB. We should point out that our fittings are only one of the data 
combinations and different combinations can certainly lead to somewhat different fitting results, on the other hand 
given the data set we use our results are robust. 

Observing inflation and dark energy dynamics are among the most important aspects of the current cosmological 
research. Currently we are not yet able to detect the dynamics on either sector and in particular the possible 
contaminations of the CMB observations on the largest scales will affect our results qualitatively^^. However our 
main concentration is to make a first detailed study in the probe of possible correlations between inflation and dark 
energy. If one starts always from the concordance ACDM cosmology one cannot achieve more subtle physics beyond 
that. And as we have shown in the cases one assume there is dynamics in only one of the sectors this can lead to 
some bias in our probe of the dynamics. If such result can be confirmed at great significance by future observations, 
it will also shed light on the study of inherent relationship between inflation and dark energy. 

Recently the WMAP three year data (WMAP3) has been released [Ztj UtI liil H9I l5fj| . with significant improvements 
on the quality WMAP3 team claims a nontrivial preference of the red-tilted power spectrum, which is different 
from the first year results and we can hopefully get some different fitting results. However our qualitative results 
depend strongly on the largest scales of CMB and as implied recently in Ref. the pixel-based likelihood of WMAP3 
on the largest scales seems to take too large a weight on the determination of cosmological parameters and one might 
take the risk of getting biased, in the present paper we do not include the new WMAP3 data for a reanalysis and 
leave it for future investigations. 
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